Whole genome sequencing analysis of alpaca suggests TRPV3 as a candidate gene for the suri phenotype.

BACKGROUND: Alpaca is a domestic South American camelid probably arising from the domestication of two wild camelids, the vicugna and the guanaco. Two phenotypes are described for alpaca, known as huacaya and suri. Huacaya fleece is characterized by compact, soft, and highly crimped fibers, while suri fleece is longer, straight, less crimped, and lustrous. The gene variants determining these phenotypes are still unknown, although previous studies suggested a dominant inheritance of the suri. Based on that, the aim of this study was the identification of the gene variants determining alpaca coat phenotypes through whole genome sequencing (WGS) analysis. RESULTS: The sample used includes two test-cross alpaca families, suri × huacaya, which produced two offspring, one with the suri phenotype and one with the huacaya phenotype. The analyzed sample was expanded through the addition of WGS data from six vicugnas and six guanacos; this because we assumed the absence of the gene variants linked to the suri phenotype in these wild species. The analysis of gene variant segregation with the suri phenotype, coupled with the filtering of gene variants present in the wild species, disclosed the presence in all the suri samples of a premature termination codon (PTC) in TRPV3 (transient receptor potential cation channel subfamily V member 3), a gene known to be involved in hair growth and cycling, thermal sensation, cold tolerance and adaptation in several species. Mutations in TRPV3 were previously associated with the alteration of hair structure leading to an impaired formation of the hair canal and the hair shaft in mouse. This PTC in TRPV3, due to a G > T substitution (p.Glu475*), results in a loss of 290 amino acids from the canonical translated protein, plausibly leading to a physiological dysfunction. CONCLUSION: The present results suggest that the suri phenotype may arise from a TRPV3 gene variant which may explain some of the suri features such as its longer hair fibre with lower number of cuticular scales compared to huacaya.

Clostridium lamae sp. nov., a novel bacterium isolated from the fresh feces of alpaca.

A novel Gram-positive, anaerobic, nonspore-forming, rod-shaped bacterium, designated strain NGMCC 1.200840 T, was isolated from the alpacas fresh feces. The taxonomic position of the novel strain was determined using a polyphasic approach. Phylogenetic analysis based on 16S rRNA gene sequences revealed strain NGMCC 1.200840 T was a member of the genus Clostridium and closely related to Clostridium tertium DSM 2485 T (98.16% sequence similarity). Between strains NGMCC 1.200840 T and C. tertium DSM 2485 T, the average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) were 79.91% and 23.50%, respectively. Genomic DNA G + C content is 28.44 mol%. The strain can utilise D-glucose, D-mannitol, D-lactose, D-saccharose, D-maltose, D-xylose, L-arabinose, D-cellobiose, D-mannose, D-melezitose, D-raffinose, D-sorbitol, L-rhamnose, D-trehalose, D-galactose and Arbutin to produce acid. The optimal growth pH was 7, the temperature was 37 °C, and the salt concentration was 0-0.5% (w/v). The major cellular fatty acids (> 10%) included iso-C15:0, anteiso-C15:0 and iso-C17:0 3-OH. The polar lipids consisted of diphosphatidylglycerol, phosphatidylglycerol, phosphatidylethanolamine, three unidentified phospholipids and two unidentified aminolipids. Based on phenotypic, phylogenetic and chemotaxonomic characteristics, NGMCC 1.200840 T represents a novel species within the genus Clostridium, for which the named Clostridium lamae sp. nov. is proposed. The type strain is NGMCC 1.200840 T (= CGMCC 1.18014 T = JCM 35704 T).

Identification of genomic regions associated with differences in fleece type in Huacaya and Suri alpacas (Vicugna pacos).

The difference in fleece type is the distinguishing trait between the two types of alpacas (Vicugna pacos), Huacaya and Suri. The Suri fleece type has been found to be inherited dominantly over the Huacaya type, resulting in offspring with the Suri phenotype. The aim of our study was to map genomic regions associated with the two different fleece types. In this study, 91 alpacas (54 Huacayas and 37 Suris) from Germany and Switzerland were genotyped using the 76k alpaca SNP array. Only 59k chromosome-localised markers map to the alpaca reference assembly VicPac3.1, and after quality control 49 866 SNPs, were retained for population structure assessment and to conduct a genome-wide association study. Both principal component and neighbour-joining tree analysis showed that the two fleece-type cohorts overlapped rather than forming two distinct clusters. Genome-wide significantly associated markers were observed in the scaffold region of chromosome 16 (NW_021964192.1), which contains a cluster of keratin genes. A haplotype predominantly found in Suri alpacas has been identified which supports dominant inheritance. Variant filtering of nine whole-genome sequenced alpacas from both fleece types in the critical interval of 0.4 Mb did not reveal perfect segregation of either fleece type for specific variants. To our knowledge, this is the first study to use the recently developed species-specific SNP array to identify genomic regions associated with differences in fleece type in alpacas. There are still some limitations, such as the preliminary status of the reference assembly and the incomplete annotation of the alpaca genome.

Genetic variability among and within domestic Old and New World camels at the α-lactalbumin gene (LALBA) reveals new alleles and polymorphisms responsible for differential expression.

α-Lactalbumin (α-LA), which is encoded by the LALBA gene, is a major whey protein that binds to Ca2+ and facilitates lactose synthesis as a regulatory subunit of the synthase enzyme complex. In addition, it has been shown to play central roles in immune modulation, cell-growth regulation, and antimicrobial activity. In this study, a multitechnical approach was used to fully characterize the LALBA gene and its variants in both coding and regulatory regions for domestic camelids (dromedary, Bactrian camel, alpaca, and llama). The gene analysis revealed a conserved structure among the camelids, but a slight difference in size (2,012 bp on average) due to intronic variations. Promoters were characterized for the transcription factor binding sites (11 found in total). Intraspecies sequence comparison showed 36 SNPs in total (2 in the dromedary, none in the Bactrian camel, 22 in the alpaca, and 12 in the llama), whereas interspecies comparison showed 86 additional polymorphic sites. Eight SNPs were identified as trans-specific polymorphisms, and 2 of them (g.112A>G and g.1229A>G) were particularly interesting in the New World camels. The first creates a new binding site for transcription factor SP1. An enhancing effect of the g.112G variant on the expression was demonstrated by 3 independent pGL3 gene reporter assays. The latter is responsible for the p.78Ile>Val AA replacement and represents novel allelic variants (named LALBA A and B). A link to protein variants has been established by isoelectric focusing (IEF), and bioinformatics analysis revealed that carriers of valine (g.1229G) have a higher glycosylation rate. Genotyping methods based on restriction fragment length polymorphism (PCR-RFLP) were set up for both SNPs. Overall, adenine was more frequent (0.54 and 0.76) at both loci. Four haplotypes were found, and the AA and GA were the most common with a frequency of 0.403 and 0.365, respectively. Conversely, a putative biological gain characterizes the haplotype GG. Therefore, opportunities for rapid directional selection can be realized if this haplotype is associated with favorable milk protein properties. This study adds knowledge at the gene and protein level for α-LA (LALBA) in camelids and importantly contributes to a relatively unexplored research area in these species.

MC1R Gene Variants and Their Relationship with Coat Color in South American Camelids.

In domestic camelids, fleece color is an essential characteristic because it defines the direction of production. Variants were determined in the MC1R gene that showed a relationship with coat color in alpacas and llamas at the level of the coding region. This report sequenced the MC1R gene from 290 alpacas (142 white, 84 black, 50 brown, and 14 light fawn), five brown llamas, nine vicuñas, and three guanacos to analyze the association between coat color and the MC1R gene among South American camelids. A total of nineteen polymorphisms were identified. Seven polymorphisms were significant; three of them were of nonsynonymous type (c.82A > G, c.376G > A, and c.901C > T), two were of synonymous type (c.126 T > C and c.933G > A), one was in the promoter region (-42C > G), and one was in the 3' UTR (+5T > C). More polymorphisms were found in domestic camelids than in wild camelids. Besides polymorphism, the association of polymorphisms might cause white and dark pigmentation in the fleece of South American camelids. In addition, the MC1R protein would answer the pigmentation in alpacas.

Genome-wide scan for runs of homozygosity in South American Camelids.

BACKGROUND: Alpaca (Vicugna pacos), llama (Lama glama), vicugna (Vicugna vicugna) and guanaco (Lama guanicoe), are the camelid species distributed over the Andean high-altitude grasslands, the Altiplano, and the Patagonian arid steppes. Despite the wide interest on these animals, most of the loci under selection are still unknown. Using whole-genome sequencing (WGS) data we investigated the occurrence and the distribution of Runs Of Homozygosity (ROHs) across the South American Camelids (SACs) genome to identify the genetic relationship between the four species and the potential signatures of selection. RESULTS: A total of 37 WGS samples covering the four species was included in the final analysis. The multi-dimensional scaling approach showed a clear separation between the four species; however, admixture analysis suggested a strong genetic introgression from vicugna and llama to alpaca. Conversely, very low genetic admixture of the guanaco with the other SACs was found. The four species did not show significant differences in the number, length of ROHs (100-500 kb) and genomic inbreeding values. Longer ROHs (> 500 kb) were found almost exclusively in alpaca. Seven overlapping ROHs were shared by alpacas, encompassing nine loci (FGF5, LOC107034918, PRDM8, ANTXR2, LOC102534792, BSN, LOC116284892, DAG1 and RIC8B) while nine overlapping ROHs were found in llama with twenty-five loci annotated (ERC2, FZD9, BAZ1B, BCL7B, LOC116284208, TBL2, MLXIPL, PHF20, TRNAD-AUC, LOC116284365, RBM39, ARFGEF2, DCAF5, EXD2, HSPB11, LRRC42, LDLRAD1, TMEM59, LOC107033213, TCEANC2, LOC102545169, LOC116278408, SMIM15, NDUFAF2 and RCOR1). Four overlapping ROHs, with three annotated loci (DLG1, KAT6B and PDE4D) and three overlapping ROHs, with seven annotated genes (ATP6V1E1, BCL2L13, LOC116276952, BID, KAT6B, LOC116282667 and LOC107034552), were detected for vicugna and guanaco, respectively. CONCLUSIONS: The signatures of selection revealed genomic areas potentially selected for production traits as well as for natural adaptation to harsh environment. Alpaca and llama hint a selection driven by environment as well as by farming purpose while vicugna and guanaco showed selection signals for adaptation to harsh environment. Interesting, signatures of selection on KAT6B gene were identified for both vicugna and guanaco, suggesting a positive effect on wild populations fitness. Such information may be of interest to further ecological and animal production studies.

Genetic parameters for fleece uniformity in alpacas.

Fiber diameter is the main selection objective and criterion in alpaca breeding programs, but it can vary across anatomic regions of the animal. As fiber diameter is usually registered from a unique sample from the mid side of the body, fiber diameter variability within fleece is never addressed and phenotypic and genetic differences may exist for fleece uniformity in alpaca populations. The objective of this work was to estimate the genetic parameters of fleece uniformity in an alpaca population. Fiber diameters measured in three different locations were used as repeated records of the same animal and studied for fitting a model that considers heterogeneous the residual variance of the model. Also, the logarithm of the standard deviation of the three measures was used as a measure of the fleece variability. Estimate of the additive genetic variance of the environmental variability was 0.43±0.14, enough high to suggest the existence of wide room to select for fleece uniformity. Genetic correlation of the trait with its environmental variability was 0.76±0.13 showing that fleece uniformity will be indirectly selected when aiming to reduce the fiber diameter. In the light of these parameters, and due to the cost of registering and the cost of opportunity, it looks no worthy to include uniformity as a selection criterion in alpaca breeding programs.

The quality of alpaca fiber is mainly assessed by a low fiber diameter. However, the fiber diameter can greatly vary along the different body locations of the animal, the industry demands not only the fineness of the fiber but also the fleece uniformity. This work studied the genetic parameters related to fleece uniformity by analyzing the diameter of three samples from different body locations (mid side, shoulder, and thigh) under two different models of analysis. The results showed variability between sampling locations and the existence of important genetic variability susceptible to being used in alpaca selection. Nevertheless, selection based on a single measurement could be used as the high correlations between locations and between the trait and its variability, saving the cost of sample analyses, being the fleece uniformity indirectly selected.

Characterization of heavy-chain antibody gene repertoires in Bactrian camels.

Camelids are the only mammals that can produce functional heavy-chain antibodies (HCAbs). Although HCAbs were discovered over 30 years ago, the antibody gene repertoire of Bactrian camels remains largely underexplored. To characterize the diversity of variable genes of HCAbs (VHHs), germline and rearranged VHH repertoires are constructed. Phylogenetics analysis shows that all camelid VHH genes are derived from a common ancestor and the nucleotide diversity of VHHs is similar across all camelid species. While species-specific hallmark sites are identified, the non-canonical cysteines specific to VHHs are distinct in Bactrian camels and dromedaries compared with alpacas. Though low divergence at the germline repertoire between wild and domestic Bactrian camels, higher expression of VHHs is observed in some wild Bactrian camels than that of domestic ones. This study not only adds our understanding of VHH repertoire diversity across camelids, but also provides useful resources for HCAb engineering.

Complete Genome Sequencing Reveals Unusual Equine Rotavirus A of Bat Origin from India.

Rotaviruses are the most common viral agents associated with foal diarrhea. Between 2014 and 2017, the annual prevalence of rotavirus in diarrheic foals ranged between 18 and 28% in Haryana (India). Whole-genome sequencing of two equine rotavirus A (ERVA) isolates (RVA/Horse-wt/IND/ERV4/2017 and RVA/Horse-wt/IND/ERV6/2017) was carried out to determine the genotypic constellations (GCs) of ERVAs. The GCs of both the isolates were G3-P[3]-I8-R3-C3-M3-A9-N3-T3-E3-H6, a unique combination reported for ERVAs so far. Both the isolates carried VP6 of genotype I8, previously unreported from equines. Upon comparison with RVAs of other species, the GC of both isolates was identical to that of a bat rotavirus strain, MSLH14, isolated from China in 2012. The nucleotide sequences of the genes encoding VP3, NSP2, and NSP3 shared >95.81% identity with bat RVA strains isolated from Africa (Gabon). The genes encoding VP1, VP2, VP7, NSP1, and NSP4 shared 94.82% to 97.12% nucleotide identities with the human strains which have zoonotic links to bats (RCH272 and MS2015-1-0001). The VP6 genes of both strains were distinct and had the highest similarity of only 87.08% with that of CMH222, a human strain of bat origin. The phylogenetic analysis and lineage studies revealed that VP7 of both isolates clustered in a new lineage (lineage X) of the G3 genotype with bat, human, and alpaca strains. Similarly, VP4 clustered in a distinct P[3] lineage. These unusual findings highlight the terra incognita of the genomic diversity of equine rotaviruses and support the need for the surveillance of RVAs in animals and humans with a "one health" approach. IMPORTANCE Rotaviruses are globally prevalent diarrheal pathogens in young animals including foals, piglets, calves, goats, sheep, cats, and dogs along with humans. The genome of rotaviruses consists of 11 segments, which enables them to undergo reshuffling by reassortment of segments from multiple species during mixed infections. In this study, the prevalence of equine rotaviruses was 32.11% in organized equine farms of North India. The complete genome analysis of two ERVA isolates revealed an unusual genomic constellation, which was previously reported only in a bat RVA strain. A segment-wise phylogenetic analysis revealed that most segments of both isolates were highly similar either to bat or to bat-like human rotaviruses. The occurrence of unusual bat-like rotaviruses in equines emphasizes the need of extensive surveillance of complete genomes of both animal and human rotaviruses with a "one health" approach.

An alpaca single-domain antibody (VHH) phage display library constructed by CDR shuffling provided high-affinity VHHs against desired protein antigens.

Antigen-combining sites of the camelid heavy-chain antibody variable domain (VHH) are constructed by three complementarity-determining regions (CDR1, CDR2 and CDR3). We prepared cDNA using mRNA extracted from peripheral lymphocytes of alpacas that had been non-immunized or immunized with human serum albumin (HSA). The VHH gene fragments encoding the amino-terminal half-containing CDR1 as well as CDR2 and the carboxy-terminal half-containing CDR3 were amplified independently by PCR, and then full-length VHH gene fragments were generated by overlap extension PCR and cloned into the phagemid vector. This protocol, referred to as CDR shuffling, allowed us to construct an alpaca VHH phage display library possessing repertoires different from those naturally occurring in animals. We asked, first, whether this library was able to provide the functional VHH fragments against HSA, an immunized antigen, and obtained 29 anti-HSA VHH clones, 41% possessed KD values of lower than 10-8 M, 5 of which had KD values of 10-10 M. We also obtained VHH clones against non-immunized protein antigens such as cardiac troponin T and I, Ebola virus glycoprotein 1 and human immunoglobulin G by biopanning. We compared the amino acid sequences and affinities and found that 43% of VHHs had KD values of less than 10-8 M, although those having KD values of 10-10 M were unavailable. These results suggested that the CDR-shuffled VHH phage display library could potentially provide VHHs against non-immunized protein antigens with similar levels of affinities to those against immunized antigens.

The ASIP gene in the llama (Lama glama): Alternative transcripts, expression and relation with color phenotypes.

The Agouti gene (ASIP) is one of the most important genes for coat color determination in mammals. It has a complex structure with several promoters and alternative non-coding first exons that are transcribed into mRNAs with different 5'UTR. These mRNA isoforms regulate the temporal and spatial expression of the gene, producing diverse pigmentation patterns. Here, we studied ASIP transcriptional variants and their expression in the skin of llamas with different coat color phenotypes. We also described the ASIP locus, including promoter usage and the splicing events that originate each transcript variant. Using 5'RACE-PCR we isolated seven ASIP transcripts with alternative 5'UTR, where exons 1A, 1A', 1C, 1D, and a novel non-coding exon 1A" were identified. Additionally, new alternative spliced forms were found. The diversity of ASIP 5'UTRs is originated by a complex pattern of alternative promoter usage, multiple transcription start sites and splicing events that include exon skipping and alternative 3' splicing site selection. We found that ASIP was highly expressed in llamas with white and brown phenotypes while black animals presented very low expression. The main responsible for this difference was a fusion transcript between ASIP and NCOA6 genes, which was present in the skin of white and brown llamas but not in the black ones. The rest of ASIP transcripts presented very low expression in the skin, indicating that the main regulation point for ASIP gene expression is at the transcriptional level. Nevertheless, the characteristics of the 5'UTRs sequences suggest that alternative transcripts could be regulated differently at the protein synthesis level.

Nanobodies from camelid mice and llamas neutralize SARS-CoV-2 variants.

Since the start of the COVID-19 pandemic, SARS-CoV-2 has caused millions of deaths worldwide. Although a number of vaccines have been deployed, the continual evolution of the receptor-binding domain (RBD) of the virus has challenged their efficacy. In particular, the emerging variants B.1.1.7, B.1.351 and P.1 (first detected in the UK, South Africa and Brazil, respectively) have compromised the efficacy of sera from patients who have recovered from COVID-19 and immunotherapies that have received emergency use authorization1-3. One potential alternative to avert viral escape is the use of camelid VHHs (variable heavy chain domains of heavy chain antibody (also known as nanobodies)), which can recognize epitopes that are often inaccessible to conventional antibodies4. Here, we isolate anti-RBD nanobodies from llamas and from mice that we engineered to produce VHHs cloned from alpacas, dromedaries and Bactrian camels. We identified two groups of highly neutralizing nanobodies. Group 1 circumvents antigenic drift by recognizing an RBD region that is highly conserved in coronaviruses but rarely targeted by human antibodies. Group 2 is almost exclusively focused to the RBD-ACE2 interface and does not neutralize SARS-CoV-2 variants that carry E484K or N501Y substitutions. However, nanobodies in group 2 retain full neutralization activity against these variants when expressed as homotrimers, and-to our knowledge-rival the most potent antibodies against SARS-CoV-2 that have been produced to date. These findings suggest that multivalent nanobodies overcome SARS-CoV-2 mutations through two separate mechanisms: enhanced avidity for the ACE2-binding domain and recognition of conserved epitopes that are largely inaccessible to human antibodies. Therefore, although new SARS-CoV-2 mutants will continue to emerge, nanobodies represent promising tools to prevent COVID-19 mortality when vaccines are compromised.

Development of a 76k Alpaca (Vicugna pacos) Single Nucleotide Polymorphisms (SNPs) Microarray.

Small farm producers' sustenance depends on their alpaca herds and the production of fiber. Genetic improvement of fiber characteristics would increase their economic benefits and quality of life. The incorporation of molecular marker technology could overcome current limitations for the implementation of genetic improvement programs. Hence, the aim of this project was the generation of an alpaca single nucleotide polymorphism (SNP) microarray. A sample of 150 Huacaya alpacas from four farms, two each in Puno and Cerro de Pasco were used for SNP discovery by genotyping by sequencing (GBS). Reduced representation libraries, two per animal, were produced after DNA digestion with ApeK1 and double digestion with Pst1-Msp1. Ten alpaca genomes, sequenced at depths between 12× to 30×, and the VicPac3.1 reference genome were used for read alignments. Bioinformatics analysis discovered 76,508 SNPs included in the microarray. Candidate genes SNPs (302) for fiber quality and color are also included. The microarray SNPs cover 90.5% of the genome length with a density of about 39 ± 2.51 SNPs/Mb of DNA at an average interval of 26.45 ± 18.57 kbp. The performance was evaluated by genotyping 30 family trios and comparing them to their pedigrees, as well as comparing microarray to GBS genotypes. Concordance values of 0.93 and 0.94 for ApeK1 and Pst1-Msp1 generated SNPs were observed. Similarly, 290 fiber quality and color candidate gene SNPs were validated. Availability of this microarray will facilitate genome-wide association studies, marker-assisted selection and, in time, genomic selection.

Ancient DNA reveals the lost domestication history of South American camelids in Northern Chile and across the Andes.

The study of South American camelids and their domestication is a highly debated topic in zooarchaeology. Identifying the domestic species (alpaca and llama) in archaeological sites based solely on morphological data is challenging due to their similarity with respect to their wild ancestors. Using genetic methods also presents challenges due to the hybridization history of the domestic species, which are thought to have extensively hybridized following the Spanish conquest of South America that resulted in camelids slaughtered en masse. In this study, we generated mitochondrial genomes for 61 ancient South American camelids dated between 3,500 and 2,400 years before the present (Early Formative period) from two archaeological sites in Northern Chile (Tulán-54 and Tulán-85), as well as 66 modern camelid mitogenomes and 815 modern mitochondrial control region sequences from across South America. In addition, we performed osteometric analyses to differentiate big and small body size camelids. A comparative analysis of these data suggests that a substantial proportion of the ancient vicuña genetic variation has been lost since the Early Formative period, as it is not present in modern specimens. Moreover, we propose a domestication hypothesis that includes an ancient guanaco population that no longer exists. Finally, we find evidence that interbreeding practices were widespread during the domestication process by the early camelid herders in the Atacama during the Early Formative period and predating the Spanish conquest.

Effect of Humanizing Mutations on the Stability of the Llama Single-Domain Variable Region.

In vivo clinical applications of nanobodies (VHHs) require molecules that induce minimal immunoresponse and therefore possess sequences as similar as possible to the human VH domain. Although the relative sequence variability in llama nanobodies has been used to identify scaffolds with partially humanized signature, the transformation of the Camelidae hallmarks in the framework2 still represents a major problem. We assessed a set of mutants in silico and experimentally to elucidate what is the contribution of single residues to the VHH stability and how their combinations affect the mutant nanobody stability. We described at molecular level how the interaction among residues belonging to different structural elements enabled a model llama nanobody (C8WT, isolated from a naïve library) to be functional and maintain its stability, despite the analysis of its primary sequence would classify it as aggregation-prone. Five chimeras formed by grafting CDRs isolated from different nanobodies into C8WT scaffold were successfully expressed as soluble proteins and both tested clones preserved their antigen binding specificity. We identified a nanobody with human hallmarks that seems suitable for humanizing selected camelid VHHs by grafting heterologous CDRs in its scaffold and could serve for the preparation of a synthetic library of human-like single domains.

An 8.22 Mb Assembly and Annotation of the Alpaca (Vicugna pacos) Y Chromosome.

The unique evolutionary dynamics and complex structure make the Y chromosome the most diverse and least understood region in the mammalian genome, despite its undisputable role in sex determination, development, and male fertility. Here we present the first contig-level annotated draft assembly for the alpaca (Vicugna pacos) Y chromosome based on hybrid assembly of short- and long-read sequence data of flow-sorted Y. The latter was also used for cDNA selection providing Y-enriched testis transcriptome for annotation. The final assembly of 8.22 Mb comprised 4.5 Mb of male specific Y (MSY) and 3.7 Mb of the pseudoautosomal region. In MSY, we annotated 15 X-degenerate genes and two novel transcripts, but no transposed sequences. Two MSY genes, HSFY and RBMY, are multicopy. The pseudoautosomal boundary is located between SHROOM2 and HSFY. Comparative analysis shows that the small and cytogenetically distinct alpaca Y shares most of MSY sequences with the larger dromedary and Bactrian camel Y chromosomes. Most of alpaca X-degenerate genes are also shared with other mammalian MSYs, though WWC3Y is Y-specific only in alpaca/camels and the horse. The partial alpaca Y assembly is a starting point for further expansion and will have applications in the study of camelid populations and male biology.

Structure- and sequence-based design of synthetic single-domain antibody libraries.

Single-domain antibody fragments known as VHH have emerged in the pharmaceutical industry as useful biotherapeutics. These molecules, which are naturally produced by camelids, share the characteristics of high affinity and specificity with traditional human immunoglobulins, while consisting of only a single heavy chain. Currently, the most common method for generating VHH is via animal immunization, which can be costly and time-consuming. Here we describe the development of a synthetic VHH library for in vitro selection of single domain binders. We combine structure-based design and next-generation sequencing analysis to build a library with characteristics that closely mimic the natural repertoire. To validate the performance of our synthetic library, we isolated VHH against three model antigens (soluble mouse PD-1 ectodomain, amyloid-ß peptide, and MrgX1 GPCR) of different sizes and characteristics. We were able to isolate diverse binders targeting different epitopes with high affinity (as high as 5 nM) against all three targets. We then show that anti-mPD-1 binders have functional activity in a receptor blocking assay.

Were domestic camelids present on the prehispanic South American agricultural frontier? An ancient DNA study.

The southern boundary of prehispanic farming in South America occurs in central Mendoza Province, Argentina at approximately 34 degrees south latitude. Archaeological evidence of farming includes the recovery of macrobotanical remains of cultigens and isotopic chemistry of human bone. Since the 1990s, archaeologists have also hypothesized that the llama (Lama glama), a domesticated South American camelid, was also herded near the southern boundary of prehispanic farming. The remains of a wild congeneric camelid, the guanaco (Lama guanicoe), however, are common in archaeological sites throughout Mendoza Province. It is difficult to distinguish bones of the domestic llama from wild guanaco in terms of osteological morphology, and therefore, claims that llama were in geographic areas where guanaco were also present based on osteometric analysis alone remain equivocal. A recent study, for example, claimed that twenty-five percent of the camelid remains from the high elevation Andes site of Laguna del Diamante S4 were identified based on osteometric evidence as domestic llama, but guanaco are also a likely candidate since the two species overlap in size. We test the hypothesis that domesticated camelids occurred in prehispanic, southern Mendoza through analysis of ancient DNA. We generated whole mitochondrial genome datasets from 41 samples from southern Mendoza late Holocene archaeological sites, located between 450 and 3400 meters above sea level (masl). All camelid samples from those sites were identified as guanaco; thus, we have no evidence to support the hypothesis that the domestic llama occurred in prehispanic southern Mendoza.